Monte Carlo study of the superfiuid weight in doped antiferromagnets 
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The phase fluctuations of the condensate in doped antiferromagnets, described by a t-t'-J model and 
a suitable 1/N expansion, provide a mechanism for a Kosterlitz-Thouless (KT) type of transition to a 
superconducting state below T c . In this paper, we present a Monte Carlo study of the corresponding 
superfiuid weight D S (T) in the classical (large- iV) limit, as a function of temperature and doping. 
Consistent with generic experimental trends, D S (T) exhibits a T- linear decrease at low temperatures, 
with the magnitude of the slope D' s (0) increasing upon doping. Finite-size scaling in the underdopcd 
regime predicts values for the dimensionless ratio A = }cbT c / D s (0) of order unity, with A = 0.4435(5) 
in the half-filled-band limit, thus confirming D s (0) as the fundamental energy scale determining T c . 
Our Monte Carlo results for D s (T)/D s (0) vs k B T/D s (0), at 10% hole doping, are found to be in 
reasonable agreement with recent measurements on La2- a; Sr I Cu04, with x = 0.10, throughout the 
temperature range below the theoretical KT transition temperature T c . 

PACS numbers: 74.40. +k, 74.20.Mn, 74.25.Bt 
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I. INTRODUCTION 

The quest on the nature of the superconducting tran- 
sition observed in doped antiferromagnets, such as the 
lamellar high-T c copper-oxides- has entered a new era. 
Indeed, a consensus is emergingp on the importance of the 
zero-temperature superfiuid weight D s (0) as the funda- 
mental energy scale determining the corresponding tran- 
sition temperature T c . An early clue for this result came 
from the empirical Uemura relation,B valid for under- 
doped cuprates, which shows roughly a proportionality 
between the independently measurable quantities T c and 
D s (0), 



k B T c = AD 8 (0) , 



(1) 



where A is a dimensionless parameter of order unity. 
Given that the superfiuid weight D s (T) , at a general tem- 
perature T, measures the "phase stiffness" of the conden- 
sate, i.e., the energy cost to produce spatial variations of 
its phase, Emery and KivelsonB argued that ([!]) provides 
strong evidence for a phase-fluctuation driven supercon- 
ducting transition. Of course, in strictly two-dimensional 
systems with continuous U(l) gauge (phase) symmetry, 
melting of long-range phase coherence with increasing 
temperature, and hence loss of superfluidity of the charge 
carriers, proceeds by thermal generation and subsequent 
unbinding of vortex- antivortex pairs. The detailed de- 
scription of this transition is given by the Kosterlitz- 
Thouless (KT) theoryoEl in the context of the classical 
XY model. A characteristic feature of the KT transition 
is the discontinuous drop of DJJ 1 ) to zero at T = T c , 
as T c is approached from below B Despite the inevitable 
rounding of this discontinuity by the weak coupling be- 
tween the copper-oxide layers, Corson et alp in a re- 
markable experiment, have provided direct evidence of 
the KT nature of the superconducting transition by ob- 
serving dynamic effects of thermally generated vortices 



in the frequency-dependent conductivity of underdoped 
cuprate thin films. _ _ 

A common feature of theoretical modelsB~t2l assum- 
ing a phase- fluctuation mechanism for the superconduct- 
ing transition, is the prediction of a T-linear decrease 
in D s (T) at temperatures well beloW|-Xt~ in accordance 
with the experimental observations .I 1 HI 12 ! However, all 
the aforementioned models, do not incorporate explic- 
itly the doping dependence of D S (T) and therefore pro- 
vide no framework for a proper treatment of the high- 
Tc cuprate superconductors as doped antiferromagnets. 
This is a major drawback since the proximity of the 
doped materials to the Mott insulating antiferromag- 
netic state at half- filling (n e = 1) leadsa to a vanishing 
D s (0) and T c , as the hole concentration (1 — n e ) tends 
to zero. The importance of the antiferromagnetic fluc- 
tuations, even well within the underdoped regime, has 
been recently confirmed in a dramatic way.pSpecifically, 
the elastic neutron-scattering measurementscj in pristine 
La2_ a; Sr a ;Cu04, with x = 0.10,0.12, have established the 
presence of static long-range antiferromagnetic order in 
the superconducting ground state, suggested by earlier 
work.EJ This conclusion is supported also by the sub- 
sequent observation of conventional two-magnon Raman 
scatteringlij in the same material. 

The above discussion underlines the need for models of 
the high-T c cuprate superconductors that can account, 
on the same footing, for the intimately coupled phase 
and spin fluctuations of the condensate. Such a model 
has been suggested some time agoEJj and shown, by one 
of the authors, to display flux quantization and a finite 
superfiuid weight, i.e., superconductivity, in its ground 
state.liZI Our model consists of a t-t'-J Hamiltonian and 
a suitable 1/N expansion that provide a framework for 
the study of: (i) the ground-state properties, using stan- 
dard "spin-wave" techniques which allow the incorpora- 
tion of leading quantum-fluctuation effects, and (ii) the 
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finite-temperature properties, using an associated clas- 
sical (large- N limit) energy functional and correspond- 
ing partition function in terms of which important phys- 
ical quantities can be readily expressed and calculated by 
Monte Caxln-aimulation. Point (i) was the subject of ear- 
lier worksEJliJ and the results, especially for the optical 
properties, were found to provide support for our effective 
model when compared with experiment. In the present 
paper we consider point (ii), focusing on the study of 
the superfluid weight as a function of temperature and 
doping, a subject of current experimental and theoretical 
interest. 

In Sec. we give a brief description of our effective 
t-t'-J Hamiltonian and the associated classical energy 
functional an d pa rtition function emerging in the large-iV 
limit. In Sec. Ill we derive an explicit expression for the 
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superfluid weight D s (T) which we study by Monte Caria 
simulation using the standard Metropolis algorithmic 
We present results for the shape of the scaled curve, 
D s (T)/D s (0) vs k B T/D s (0), the value of the dimen- 
sionless ratios fcsT c /Z) s (0) and D s (T c )/kBT c , as well as 
the value of the zero-temperature slope D' s (0), clarify- 
ing their doping dependence in the regime of interest, 
i.e., close to half-filling. The KT nature of the transi- 
tion to the superconducting state is supported by finite- 
size scaling analysis using the relevant Weber-Minnhagen 
scaling formulaJ13 In particular, our numerical results for 
the temperature dependence of D S (T), at 10% hole dop. 
ing, compare reasonably well with recent measurementsEBI 
in underdoped La2_ a; Sr a ,Cu04, for the corresponding 
concentration value x — 0.10. Furthermore, the zero- 
temperature slope of the superfluid weight is predicted 
to approach, close to half-filling (n e — > 1), the universal 
value: D' S (Q) = -^kaJ2, consistent with the available ex- 
perimental data.EH The latter limiting value is shown 
to be a hallmark of the strong antiferromagnetic corre- 
lations, present in this temperature and doping regime. 
Our concluding remarks are summarized in Sec. IV. 



t if i, j are nearest neighbors 
— t' if i,j are next nearest neighbors 
otherwise . 



(3) 



The conventions in (||) incorporate opposite signs forj-t 
and t' as it is appropriate for the hole-doped cuprates.E2l 
In Ref. [l6] we generalized the local constraint associated 
with (||) to Xi° + xf M = N, where N is an arbitrary in- 
teger, and considered the commutation properties of the 
X ab operators to be those of the generators of the U(3) 
algebra. A Holstein-Primakoff realization for the latter 
algebra in terms of hard-core bosons resolves then explic- 
itly the local constraint and can be used to develop a per- 
turbation theory based on the l/N expansion, restoring 
the physical value N = 1 at the end of the calculation. 

In the presence of an external magnetic flux thread- 
ing the two-dimensional lattice in an Aharonov-Bohm 
torus geometry, the hopping matrix elements t^ arc mod- 
ified by the well known Peierls phase factor and should 
be substituted in (0) according to 



tij tijC 



with Aij = — — — 



(Ri — Rj) ■ e x . (4) 



Here R^ is the position vector for lattice site i, e x is 
the unit vector along the ir-axis encircling the flux lines 
and $o = 2irhclq is the so-called flux quantum. As ar- 
gued in Ref. 17, in the context of the present effective 
model, whereby carriers are treated as hard-core bosons, 
the charge q entering $ should be set equal to q = 2e, 
where e is the electronic charge. 

In the large- A*" limit "condensation" occurs, i.e., the 
Bose operators become classical commuting fields. Con- 
sidering only uniform density configurations, the corre- 
sponding classicalj-epergy functional resulting from (|^)- 
H) takes the formEZI 



H($) = -n e (l-n e )^*i 



II. EFFECTIVE MODEL 

The effective model under consideration is described 
by a t-t'-J Hamiltonian expressed in terms of Hubbard 



operators X 



ah 



\a)(b\ as 



H 
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^xT-xTxT)^ (2) 



where the index corresponds to a hole, the Greek indices 
fi,i>, . . . assume two distinct values, for a spin-up and a 
spin-down electron, and the summation convention is in- 
voked. Here J is the antiferromagnetic spin-exchange in- 
teraction between nearest-neighbor sites on a square 
lattice endorsed with periodic boundary conditions and 
a total number of sites A = A x x A y , where A^ = A y . 
For the hopping matrix elements Uj we assume 



cos — cos — cos ( 



- sin — sin — cos A;* 
2 2 1 3 



i/>i ~ ipj - §i + 1 



ipi - ipj +4>i- 4>j 



(5) 



where n e is the average electronic density, the angles 
(9i,<f>i) determine the local spin direction, and the re- 
maining parameter ipi determines the local phase of the 
condensate. The above functional form makes apparent 
the coupling between the phase and spin variables of the 
condensate through the kinetic energy term, proportional 

tO tij . 

The description of the finite-temperature classical the- 
ory is now completed using the energy functional ([5]) to 
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construct the corresponding partition function Z(&) and 
free energy per lattice site F($>), 



Z(#) 



= e -/5A_F($) = 




(0) 



where = l/(fcsT) and the integrations at each lattice 
site i extend over the intervals: < 9i < tt, < (pi < 2ir, 
and < ipi < 47r. Invoking standard thermodynamic 
identities, important physical quantities can be readily 
expressed in terms of the partition function (^|) and stud- 
ied by Monte Carlo simulation. In the following we fo- 
cus our discussion on the study of the superfluid weight 
D S (T). 



III. SUPERFLUID WEIGHT 

At a finite temperature T, the superfluid weight (or 
helicity modulus) D s (T) is given by the curvature of the 
infinite lattice limit of the free energy AF($) at <f> = 



m 2 


'd 2 F(<S>)~ 




<9$ 2 



(7) 



Physically, D s (T) determines the ratio of the density of 
the superfluid charge carriers to their mass and hence 
can be related to the experimentally measurable in-plane 
magnetic (London) penetration depth, as noted later on 
in this section. Carrying out the second derivative with 
respect to $ in (|^) we have more explicitly that 



D S (T) = n e (l 



R-7 



cos — cos — cos 
2 2 



sin — sin — cos 
2 2 



ipi -ipj ~<Pi+- 



ipi ~ ipj + ■ 



Mi 



cos — cos — sin 
2 2 



+ sin — sin — sin 
2 2 




ipi - ipj + <pj 



(8) 

z — 4 being the coordination number of the square lat- 
tice. As shown in Ref. [l6], close to half-filling (n e < 1) 
and for a sufficiently large t' , the ground state of (|), in 
the absence of magnetic flux (<I> = 0), is described by a 



planar spin configuration (6i — ir/2) in which the local 
twist angles and phases exhibit long-range order accord- 
ing to: (pi = Q -Rj, ipi = Q' ■ Ri, where Q = (ir, tt) is the 
usual spin-modulating antifcrromagnetic wavevector and 
Q' = (7r, — 7r) is a phase-modulating wavevector. The 
zero-temperature value of the superfluid weight follows 
then easily from (||) as 



D s (0) = 4t'n e {l - n e ) 



(9) 



For the typical two-dimensional model with continuous 
symmetry under consideration, we expect that at low 
but finite temperatures, the long-range order will be de- 
stroyed by the proliferation of excited Goldstone modes, 
leading to a T-linear decrease of D S (T). At higher tem- 
peratures we expect that the thermal generation and 
subsequent unbinding of vortex-antivortex pairs will lead 
eventually to a discontinuous drop of D S (T) to zero, at a 
critical point T — T c , in a KT type of transition. In or- 
der to study numerically D a (T) in the whole temperature 
range and affirm the aforementioned physical picture, we 
performed a Monte Ca*lo simulation using the standard 
Metropolis algorithmLa Our calculations were carried 
out on small lattices, with typical sizes A x = 8, 16, 32, 64. 
For a given temperature we performed of the order of 10 4 
thermalization steps and of the order of 10 6 measure- 
ments. We considered values for the dimensionless ratios 
e = t'/t = 0.45 and t/J — 1.0, which are thought to be 
relevant for the copper-oxide layers, and restricted our 
study to the underdoped regime, i.e., to small (1 — n e ) 
values up to 10% hole doping. The latter restriction 
is dictated by the fact that models of the t-t'-J kind, 
being rather simple extensions of the Mott-Heisenberg 
antifcrromagnetic insulator, cannot properly account for 
the nontrivial evolution of the electronic structure of the 
cuprates that occurs at higher doping values, namely, the 
closing of the pseudogapO in the optimally doped and 
overdoped (Fermi liquid) regime. 

Typical Monte Carlo results for the superfluid weight 
vs temperature are shown in Fig. |l](a) , for (1— n e ) = 0.01, 
and Fig. |l|(b), for (1 — n e ) ~ 0.10. At low temperatures 
the superfluid weight has a weak finite-size dependence 
and displays the expected T-linear decrease. In partic- 
ular, for T — > and n e — » 1, we have established the 
asymptotic form 



Ds(T) 
D.(0) 



= 1- 



k B T 
2D S (0) 



for T^0 and 



1. (10) 



The analytic expression (10) is shown in Fig. [l]by a dot- 
ted line. Evidently, this asymptotic line is approached 
very closely from below by the low temperature numer- 
ical data already in the case of the 1% hole doping; see 
Fig. |l|(a). From @ follows that the zero-temperature 
slope of the superfluid weight approaches, close to half- 
filling, the parameter-free universal value: 



D' s (0) = -kg/2 



-0.043 meVIC 



for n e — > 1 . (11) 



3 



The upper limiting value ( |11| ) is a rather stringent pre- 
diction of our theory and seems, indeedr-to be consistent 
with the available experimental datafcHtJ in the high-T c 
cuprate superconductors. A comparison of Fig. 0(a) with 
Fig. |l|(b) reveals an increase in the magnitude of the 
slope D' s (0\-U|ppn doping, a trend also consistent with 
experiment Ol 
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FIG. 1. Superfluid weight D S (T) vs temperature, for var- 
ious lattice sizes, e = 0.45, t/J = 1.0, and (a) 1 — n e = 0.01 
[estimated k B T c /D a (0) = 0.4346(9)], (b) 1 - n e = 0.10 [esti- 
mated k B T c /D s (0) = 0.3502(8)]. Monte Carlo points above 
the corresponding estimated KT transition temperature T c 
are nonzero due to finite-size effects. Error bars are included 
but in most cases are smaller than symbol size. 



case, the vanishing overlap between the opposite sublat- 
tice spin states leaves the direct hopping t' between next- 
nearest-neighbor (n.n.n) sites as the only relevant process 
of charge transport. One can then easily show that the 
expressions for the energy functional ^ and superfluid 
weight (||) reduce to those of a classical XY model for the 
ipi variables, but with only a n.n.n interaction /, where 
/ = D s (0)/2. The reduction of the structure of the phase 
fluctuations, close to half-filling, to that pf the, n.n.n XY 
model and not to the commonly assumecHL3£j structure 
of the nearest-neighbor (n.n) XY model, is an important 
prediction and a direct consequence, in the context of 
our theory, of the presence of strong antifcrromagnetic 
correlations in this regime. Numerically, the validity of 
our argument becomes apparent in Fig. |2j displaying al- 
most coinciding D S (T) Monte Carlo data for the n.n.n 
XY model (opaque diamonds) and the t-t'-J model with 
a very small hole concentration (1 — n e ) = 0.01 (filled di- 
amonds). Therefore, we may exploit the detailed results 
of Appendix |A| for the n.n.n XY model, see Eq. (|A5|), to 
conclude the asymptotic form ( |Tc| ) and hence the limit- 
ing value (pT|). The latter value, D' S (Q) = — fcn/2, being 
twice that of the n.n XY model (see Appendix |A|) , serves 
as a distinct hallmark of the sublattice structure of the 
strong antifcrromagnetic correlations in the limit n e — > 1 . 
Our observations here affirm also, by analogy to the well- 
known physics of the XY model, the presence of a KT 
transition for the superfluid weight of the t-t'-J model, 
when n e — > 1, and allow to transcribe relevant results 
for the former model, established in Appendix [a], to the 
latter, e.g., 



Ds(T c ) 4 

= — for n e — > 1 , 

k B T c 7T 

k B T c /D s {0) = 0.4435(5) for 



(12) 



(13) 



In Fig. and Fig. | we present (4/n)k B T/D s (0) by a 
short dashed line. According to (p^), in the limit n e — > 1 
the latter line should intersect the corresponding Monte 
Carlo data curve of D s (T) /D a (0) vsk B T/D s (0) precisely 
at the fcsT c /D s (0)value given by (fi"3|); see opaque (filled) 
diamonds in Fig. H. 



We emphasize that the asymptotic form (|1C]) is a phys- 
ical consequence of the fact that close to half-filling the 
antiferromagnetic exchange energy, Jrig, overwhelms the 
hole kinetic energy, ti£n e {\ — n e ), and in particular the 
term D s (0) given by (O). Hence, in the limit n e — > 1 and 
in the relevant temperature range k B T < D s (0) Jn 2 e 
(so that T — > 0), the thermal average (||) may be simpli- 
fied by freezing the spin variables (6i,(f>i) to their zero- 
temperature antiferromagnetic configuration and allow- 
ing fluctuations only in the phase variables ipi. In this 
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k B T/D s (0) 

FIG. 2. Superfluid weight D S (T) vs temperature. Experi- 
mental data (triangles) on I^-iSr^CuGi, with x = 0.1ft, are 
extracted from the measurements of Panagopoulos et alS3 us- 
ing Eq. (^). Corresponding Monte Carlo results (circles) are 
calculated for lattice size 64 x 64 and e = 0.45, t/J — 1.0, 
1 — n e = 0.10. Also included are results for the same lattice 
size and e = 0.45, t/J= 1.0, 1 — n e — 0.01 (filled diamonds), 
which should be compared with the results for the n.n.n XY 
model (opaque diamonds). The latter model corresponds to 
the limiting form of the present t-t'-J model, when n e — > 1, 
as discussed in the text. 

It should be noted that in the actual simulations of 
the t-t'-J model we can only use finite, though possi- 
bly small, hole concentrations which inevitably lead to 
deviations from the limiting values (|l^-(p"3"|). In order 
to obtain rather accurately the corresponding transition 
temperature T c we have used-the finite-size scaling anal- 
ysis of Weber and MinnhagenEj, which is appropriate for 
KT type of transitions. In this analysis one measures the 
chi-square values \ 2 (T) of the fitting of the Monte Carlo 
data for the superfluid weight, at each given tempera- 
ture and for a sequence of small lattice sizes, to a cer- 
tain scaling formula, derived by the latter authors from 
the Kosterlitz renormalization group equations.c3 Specif- 
ically, one assumes at each temperature T the following 
A^-dependence of the superfluid weight D S (T,A X ), 

where Ao(T) is some characteristic length of the order of 
the lattice constant and has no singularity at T = T c . 
The logarithmic lattice size dependence involved in ( |l4| ) 
is characteristic of the presence of vortices in the KT 
transition. Strictly speaking (|lj) is correct only at the 
critical point T = T c . Given now a value for R(T C ), the 
critical temperature can be determined by two alterna- 
tive procedures: (a) we fix R(T) to be i?(7U and use 
Ao(T) as the only adjustable parameter in ( jl4| ) to mea- 
sure the chi-square values X 2 (T) of the fitting, in which 
case T c corresponds to the minimum of the x 2 (T) curve, 
or (b) we use both R(T) and Aq(T) as adjustable param- 



eters in (|1J) and determine T c from the point where the 
R(T) curve crosses the line R(T) = R(T C ). The correct 
value for R(T C ) should lead uniquely to the same value for 
T c in both procedures. The application of this finite-size 
scaling to the n.n.n XY model is summarized pictorially 
in Fig. H of Appendix |X| and justifies, as far as the t-t'- 
J model is concerned, the limiting values (|L^-(p"3|) with 
R{T C ) -» 2, for n e -> 1. 

For the finite doping value (1 — n e ) — 0.01 or (1 — n e ) = 
0.10 (and fixed e = t'/t = 0.45, t/J = 1.0), the ap- 
plication of the aforementioned finite-size scaling anal- 
ysis, using the lattice size sequence A x = 4,6,8,10,12, 
leads to: k B T c /D s (0) = 0.4346(9) with R(T C ) = 2.04, 
or k B T c /D s (0) = 0.3502(8) with R{Tr) = 2.51, respec- 
tively. To be sure, the presence in Fig. [y of nonzero D s (T) 
values above the corresponding T c , instead of a discon- 
tinuous drop to zero, is due to finite-size effects which 
grow rapidly above the estimated critical temperature; 
a typical behavior for a KT transition. Indeed, follow-. 
ing an original argument by Weber and Minnhagen,E2l 
we note that the success of the scaling formula ( |l4[ ) in 
the present model provides, ipso facto, strong evidence 
that the relevant phase transition is of the KT type. 
Our results show a modest increase of the jump ratio 
D s {T c )/(k B T c ) = 2R{T c )/ir upon doping. The latter be- 
havior seems physically similar to that -known in the lit- 
erature of the frustrated XY modelED'Ea given that, in 
the context of the t-t'-J model, doping induces a form 
of dynamic frustration for the phase variables, via their 
inevitable coupling to the fluctuating spin variables. 

Furthermore, our results show that the dimensionless 
parameter A = k B T c / D s (0), introduced in context of the 
empirical Uemura relation (Q), is not doping independent 
but decreases modestly upon doping, while away from 
half- filling it also depends on the couplings e and t/J. 
Nevertheless, for rough theoretical estimates of the KT 
transition temperature T c in terms of D s (0), in the un- 
derdoped regime, one may always use the universal lim- 
iting value A — 0.4435(5), for n e — > 1, quoted in (|l^). 
We remind that that the latter value is characteristic of 
the n.n.n XY model and equals to half the correspond- 
ing value of the n.n XY modcLIsee Appendix |A|) com- 
monly employed to this endE'EZrtJ. Note that the use of 
A = 0.4435(5) in conjunction with (|l|) brings earlier theo- 
retical overestimations of the. KT. transition temperature 
for the copper-oxides layers,EZlE3 derived with A w 0.9, 
down to more reasonable values. In all cases, the present 
analysis confirms D s (0) as the fundamental energy scale 
determining T c in the underdoped regime. 

Having discussed the generic trends of the superfluid 
weight as a function of temperature and doping, in the 
underdoped regime, it is instructive to provide a more 
direct comparison of our theory with experiment. To 
this end we note that in the lamellar high-T c supercon- 
ductors, the experimental value for the superfluid weight 
per copper-oxide plane, JJ~^ p \t), can be extracted from 
the directly measurableOO in-plane magnetic (London) 
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penetration depth, A a (,(T), using the relation 



IS 



£)(exp)( T ) = 



{hcfd 



(15) 



where d is the average distance between planes and we 
remind that q — 2a_ Using the experimental data of 
Panagopoulos et al$3 for A a &(T) on the underdoped 
La2- 2; Sr 2 ;Gu04, with x = 0.10, and the structural 
parameters d = 6.64 A, we depict in Fig. @_by trian- 
gles the corresponding experimental values (p~q) for the 
superfluid weight vs temperature. As a result of the weak 
coupling between the copper-oxide layers, the experimen- 
tal data display no KT discontinuity but rather a con- 
tinuous drop of the superfluid weight to zero, at a spe- 
cific temperature value, that is not simply related to the 
ideal KT transition temperature of a copper-oxide mono- 
layer. Corresponding Monte Carlo results for 10% hole 
doping are depicted in Fig. § by circles and calculated 
for a 64 x 64 lattice, with e = 0.45 and t/J = 1.0. As 
noted earlier in this section, the theoretical KT tran- 
sition temperature for the latter set of parameters is 
k B T c /D s (0) = 0.3502(8), while nonzero Monte Carlo 
points above this value are due to finite-size effects. Evi- 
dently, our theoretical results (circles) in Fig. || compare 
reasonably well with the experimental data (triangles) 
throughout their common relevant temperature range, 
i.e., up to T c . 

In Fig. ||we also depict Monte Carlo results for 1% hole 
doping (filled diamonds), as well as results for the n.n.n 
XY model (opaque diamonds), thus providing the theo- 
retical lineshape of the superfluid weight vs temperature, 
in the limit n e — ► 1. Clearly, it will be very interesting to 
have measurements of the in-plane magnetic penetration 
depth on I^-xSr^CuO^ with hole concentration x as 
small as it is experimentally possible, to compare with 
the present definite theoretical prediction. 



n.n.n XY model. In order to check this definite theoreti- 
cal prediction we have suggested measurements of the in- 
plane magnetic penetration depth in very lightly doped 
cuprates. Here we should add that higher order 1/iV, 
corrections are expected to renormalize downwardsliH 
the fundamental energy scale D s (0) but, nevertheless, 
leave the lineshape of the scaled curve D S (T)/ D S (Q) vs 
ksT/ D s (0) essentially intact. 

The present study shares some common features with 
earlier worksOLj invoking a phase-fluctuation mecha- 
nism for the high-T c superconductivity. On the other 
hand, all these works including ours arep&adically dif- 
ferent from phenomenological approachescJ'U that im- 
plicate the thermally excited nodal quasiparticles in a 
e?-wave BCS superconducting state for the reduction of 
D S (T) with increasing temperature. The weak-coupling 
BCS type of approaches, however, are undermined by 
the absence of normal state quasiparticle peaka53 near 
the Brillouin zone points (0, ±7r) and (±7r, 0) where su- 
perconductivity is supposed to originate. At present it is 
still difficult to discern experimentally whether the tem- 
perature and doping dependence of D S {T) is dominated 
by phase fluctuations or by nodal BCS-like quasiparticle 
excitations. However, recent experiments in cuprate thin 
films have provided strong evidence for theinherent two- 
dimensional character of superconductivityEland for the 
KT nature of the superconducting transition!! 
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APPENDIX A: GENERALIZED XY MODEL 



IV. CONCLUDING REMARKS 

In this paper we have presented a study of the tem- 
perature and doping dependence of the superfluid weight 
D S (T) in doped antiferromagnets described by the t-t'-J 
model (||)-(||). Using Monte Carlo simulations and finite- 
size scaling analysis we have demonstrated that the phase 
fluctuations of the condensate, emerging in an appropri- 
ate classical (large- N) limit, drive superconductivity via 
a Kosterlitz-Thouless type of transition. Our theoretical 
results reproduce important generic experimental trends 
of D S (T), observed in the underdoped high-T c cuprate 
superconductors. This includes the T-linear decrease of 
D s (T) at low temperatures and the increase of the mag- 
nitude of the slope D' s (0) upon doping. 

In particular, the sublattice structure of the strong an- 
tiferromagnetic correlations in the half-filled-band limit 
was shown to dictate the lineshape of D S (T)/ D s (0) vs 
fcfiT '/ 'D s (0), for n e — » 1, to be identical to that of the 



In the main body of the paper we noted that the struc- 
ture of the phase fluctuations of the t-t'-J model un- 
der study reduces, in the half-filled-band limit, to that 
of a classical XY model with only next-nearest-neighbor 
(n.n.n) interactions. In order to clarify the properties of 
the latter model, in juxtaposition to those of the more 
conventional nearest-neighbor (n.n) XY model, we con- 
sider briefly in this appendix the following Hamiltonian 

Hxv = ~\^2 Z « cos( ^ ~ ^J") ' ( A1 ) 

assuming 

{(1 — a)I if i,j are nearest neighbors 
al if i,j are next nearest neighbors 
otherwise , 

(A2) 
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where a is a free parameter, with < a < 1, and / > 0. 
At each lattice site i the angle tpi varies in the interval: 
< ipi < 2ir. Evidently, a = corresponds to the n.n 
XY model, while a = 1 corresponds to the n.n.n XY 
model. 

The superfluid weig ht ( or hclicity modulus) for the 
generalized XY model ( Al ) reads 



D s (T)/D s (0) vs k B T/D s (0) for the n.n.n XY model co- 
incide, within numerical error, with those for the n.n XY 
model, when the ksT/ D s (0) axis is scaled by a factor of 
2. This agreement becomes better with increasing lattice 
size. 



zA\2 



0.6 - 



(A3) o" o.4 



in agreement with corresponding early results @ In view 
of (]A2|), the ground state configuration of ( |A1[ ) is simply 
given by = 0, while the zero-temperature value of the 
superfluid weight follows immediately from (|A3|) as 



#,(0) = (l+a)I 



(A4) 



Integrating the quadratic (Gaussian) fluctuations around 
the ground state configuration we obtain, after some 
lengthy algebra, the following low-temperature asymp- 
totic expansion for the superfluid weight 



for T -> 



(A5) 



D s (0) v "zD s {0) 

Here G(a) is a dimensionless geometric factor given by 



G(a) = ir 7 -- 

y ' A ^ 1 -a (1 



a(l - <S q ) 



(l-a)(l- 7q )+a(l-* q ) ' 



(A6) 



with 



7 q = \ (cos q x + cos g y ) , <S q = cos q x cos g y . (A7) 

We emphasize that G(a) is an increasing function of a 
with end-point values: G(0) = and G(l) = 1. Hence, 
from (A5) follows that the zero-temperature slope D' s (0) 
evolves monotonically from — fcs/4 to — fcs/2, as the pa- 
rameter a varies from (n.n XY model) to 1 (n.n.n XY 
model); see dash-dotted and dotted lines in Fig. ||. 

Noting now that |R; — Rj| 2 = 2, for n.n.n sites i,j, 
whereas |Ri — R,| 2 = 1, for n.n sites a cautious 



inspection of (A3) reveals that, in the thermodynamic 
limit and at each given temperature T, the superfluid 
weight of the n.n.n XY model should be twice as large 
the superfluid weight of the n.n XY model: D { s a - 1] (T) = 
2D l f~ \T). The latter property is ex plicit in the low- 
temperature analytic results (AA)-{ A5) and we have con- 
firmed its validity in the whole temperature range by 
Monte Carlo simulation using the standard Metropolis 
algorithm with the parameters quoted in Sec. HI. In- 
deed, as shown in Fig. 0, the Monte Carlo data curves of 




4- ^ 8x8 

5- v 16x16 
d-p 32x32 
3-o 64x64 

(4/jt)k B T/D s (0) 

- (2/it)k B T/D s (0) 
1-k B T/[2D s (0)] 

- 1-k B T/[4D s (0)] 



1.0 1.5 2.0 

k B T/D s (0) 

FIG. 3. Superfluid weight D S (T) vs temperature for the 
XY model with various lattice sizes and n.n.n or n.n interac- 
tion only. The crossings of the short (long) dashed line with 
the Monte Carlo data curves for the n.n.n (n.n) XY model 
provide estimates for the KT transition temperature T c of in- 
creasing accuracy, as the lattice size increases. Monte Carlo 
points above the corresponding estimated T c are nonzero due 
to finite-size effects. The dotted (dash-dotted) line depicts the 
low-temperatur e as ymptote for the n.n.n (n.n) XY model, ac- 
cording to Eq. ( |A5| ) . 

For the conventional n.n XY model it is well-known! 
that in the thermodynamic limit we have the jump ratio: 
D s (T c )/(kBT c ) = 2/tt, or equivalently in the notation of 
@, the value R(T L ) = 1. Hence, in view of the simple 



relation D { s a ~ 1] (T) = 2L>r" UJ (T), we anticipate for the 
n.n.n XY model the jump ratio: D s (T c )/(kBT c ) — 4/7T, 
or equivalently, R{T C ) — 2. In order to demonstrate nu- 
merically the latter property we have carried out a finite- 
size scaling analysis for the n.n.n XY model based on 
the scaling formula (|l4]), as described in detail Sec. III. 
The results are summarized in Fig. |i. We note that the 
minimum of the X 2 (^) curve in Fig. 3(a) occurs, within 
numerical error, at the same point were the R(T) curve 
in Fig. E|(b) crosses the line R{T) — 2. Hence the assign- 
ment R{T C ) — 2 leads, indeed, to a uniquely determined 
value for T c which from the crossing point in Fig. [|(b) 
is estimated to be k B T c /D s (0) = 0.4435(5). The suc- 
cess of the finite-size scaling analysis validates then the 
assignment R(T C ) — 2 for the n.n.n XY model. 



)(«=<>) 
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• A x =4,6,8,10,12 



0.46 



0.44 0.45 
k B T/D s (0) 

FIG. 4. Finite-size scaling around T c for the n.n.n XY 
model, according to Eq. ( jil] ) and for the lattice size sequence 
A x = 4, 6, 8, 10, 12. (a) Chi-square values (circles) of the fit- 
ting vs temperature with R(T) fixed to 2. The inset is an en- 
larged view of the curve around its minimum, (b) Coefficient 
R{T) vs temperature. The dotted vertical line indicates the 
critical temperature at which the solid line crosses the dotted 
horizontal line [R(T) = 2]. The solid line is determined by 
a linear fitting to the original data (circles). The estimated 
critical temperature is given by fesTc/ D s (0) = 0.4435(5). 

It is worth emphasizing that in terms of the dimen- 
sionless parameter A — ksTc/ D s (0), the Monte Carlo 
results of the present appendix imply the value A^ a=1 ^> — 
0.4435(5) for the n.n.n XY model and, of course, twice 
as large corresponding value for the n.n XY model, i.e., 
^(a=o) = 0.887(1). Note that the latter value for the 
n.n XY model agrees, with the original estimate of We- 
ber and Minnhager£3 derived with the same finite-size 
scaling procedure. Pictorially, our results are manifest in 
Fig. [| were the crossings of the short (long) dashed line 
with the Monte Carlo data curves for the n.n.n (n.n) XY 
model provide estimates for the value of A of increasing 
accuracy, as the lattice size increases. 
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